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Abstract 


The image quality of a large field survey telescope with a fast focal ratio of the primary mirror is high sensitivity to 
the optical elements’ misalignments and the primary mirror’s deformations. To maintain good optical performance, 
the perturbations need to be aligned, and the tomographic reconstruction of these perturbations can be derived from 
wavefront sensing with multi-field points for alignment. This work undertakes a comprehensive examination 
related to the implementation and optimization of a multi-field split curvature wavefront sensing system, including 
defocus distance, signal pre-processing, deblending of overlapped doughnuts, field-dependent correction, and 
distorted coordinate correction. We also conduct experiments to demonstrate the application and performance of a 
multi-field split curvature wavefront sensing system in Mephisto. In the context of the decentering of the secondary 
mirror, the coefficient of determination (R) attains a high value of 0.993, indicating a robust linearity between the 
coma coefficients and the decentering of the secondary mirror. A comparative analysis of the simulated and 
experimental results shows that the difference between them is less than 0.1 A. 
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1. Introduction 


Over the past two decades, large field survey telescopes have 
played a pivotal role in addressing a multitude of scientific 
objectives, such as creating comprehensive sky surveys to aid 
in the discovery of new astronomical phenomena (Wang 2018), 
investigating the nature of dark energy and dark matter by 
observing large-scale structures in the universe (Oguri et al. 
2012), studying the evolution and distribution of galaxies to 
understand their formation and development over cosmic time 
(Padilla & Strauss 2008) and exploring other critical aspects of 
astrophysics research. 

The imperative for enhanced image quality and higher time 
frequency in the ever-expanding landscape of scientific 
research propels the development of the next generation of 
sky survey telescopes. Within the latest generation of sky 
survey telescopes, there is a discernible trend toward increasing 
the aperture of the primary mirror to enhance the light- 
gathering capabilities. This necessitates the adoption of a fast 
focus ratio for the primary mirror to minimize the secondary 
mirror’s obscuration and shorten the barrel’s length. Exemplary 
instances, including J-PAS (Cenarro et al. 2018), LSST (Stubbs 
et al. 2004), Mephisto (Yuan et al. 2020), KDUST (Yuan et al. 
2012), and WFST (Lou et al. 2016), are characterized by a 
primary mirror focal ratio of less than 2. It is a very challenging 
task to maintain the high performance of this kind of telescope 


during the observation, primarily attributed to the heightened 
sensitivity of image quality to both the misalignment of optical 
elements and environmentally induced deformations of the 
primary mirror. Active optics, facilitated by measuring 
wavefronts at multiple field points, emerges as a viable 
solution for correcting the misalignments of optical elements 
and deformations in the primary mirror. 

In large field survey telescopes of the new generation, the 
measurement of wavefronts at multiple field points is 
predominantly facilitated by split curvature wavefront sensors 
(Xin et al. 2015; Chueca et al. 2012). Distinguished from the 
traditional Shack-Hartmann wavefront sensor (SHWFS), the 
split curvature wavefront sensor offers notable advantages in 
source selection owing to the large field of view afforded by 
area sensors. Moreover, SHWFS necessitates placement near 
the focal plane to avoid vignetting in a fast focal ratio system. 
This requirement imposes constraints on the dimensions of 
reimaging optics and micro-lens arrays for effective encapsula- 
tion within the CCD’s vacuum cavity. In contrast, the split 
curvature wavefront sensor presents a comparatively uncom- 
plicated configuration, facilitating integration into the CCD 
panel. 

This work undertakes a comprehensive examination related 
to the implementation and optimization of a multi-field split 
curvature wavefront sensing system, such as defocus distance, 
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Figure 1. The optical configuration of the curvature wavefront sensor. 


Table 1 
Curvature Wavefront Sensor for Sky Survey Telescopes 

Implementation DOF of the Field Defocus 
Telescope F/# of CWFS Perturbations Points Length 
VISTA 3.6 Beam splitter 5 2 1.2mm 
J-PAS 3.5 Split CWFS 7 4 1 mm 
DECam 2.7 Split CWFS 5 8 1.5 mm 
LSST 1.23 Split CWFS 24 4 1 mm 
WFST 2.48 Split CWFS 14 4 1 mm 


pre-processing of signal sources, the decoupling of signals 
from neighboring sources, field-dependent corrections, and 
distorted coordinate correction for intensity compensation. 
Subsequently, the curvature wavefront sensing system on 
Mephisto is evaluated through the experiment, and the 
simulation of active alignment procedures is tested. 


2. Principle of Multi-field Split Curvature Wavefront 
Sensing 


2.1. Configuration of the Multi-field Split Curvature 
Wavefront Sensor 


For the traditional curvature sensor, the acquisition of intra- 
focus and extra-focus images typically involves the use of a 
beam splitter (Clark et al. 2004) or physically moving the 
detector (Roddier & Roddier 1993). However, the application 
of a beam splitter strategy encounters challenges for integrating 
into the crowded focal plane, as the volume of the beam splitter 
must increase to avoid the vignetting for the fast beams. 
Alternatively, the method employing a distorted diffraction 
grating (Fengjie et al. 2007) or a non-coaxial Fibonacci-sieve 
filter (Zhang et al. 2019) can capture both intra-focus and extra- 
focus images in a single-shot. Nevertheless, these elements 
pose practical design challenges in dealing with sources in 
random fields. Notably, the split curvature wavefront sensor 
(CWFS) has gained prevalence in sky survey telescopes, as 


detailed in Table 1. The intra- and extra-focal sources of the 
split CWFS come from two different parts of the field of view, 
each captured by separate CCDs: one positioned in front of the 
focus and the other situated behind it. 

Multi-field split curvature wavefront sensors are typically 
positioned at the marginal fields of the focal plane. This 
strategy not only preserves fields characterized by high image 
quality for science but also reduces the condition number of the 
sensitive matrix and enhances the stability of tomographic 
reconstruction for the perturbations. The scalability of the 
sensitive matrix’s dimension is achieved by augmenting the 
number of measured field points, offering a more robust 
solution to the perturbation solving process. Drawing on a 
simulation conducted by Liang et al. (Liang et al. 2012) for the 
Large Synoptic Survey Telescope (LSST), a minimum of three 
curvature wavefront sensors is necessary to ensure the sensitive 
matrix attains full rank for the calibration of the 24 perturbation 
parameters (10 rigid body motions and 14 mirror deforma- 
tions). Furthermore, employing wavefront sensors at four field 
points proves effective in reducing the condition number of the 
sensitive matrix. 


2.2. Wavefront Estimating Algorithms 


In the context of the paraxial approximation, the intensity 
transport equation (ITE) quantitatively describes the intricate 
relationship between the wavefront of the beam and intensity 
variations in light intensity along the optical axis (Streibl 1984). 
As the optical configuration of the curvature wavefront sensor 
depicted in Figure 1, under the assumption of a uniformly 
illuminated pupil with intensity denoted as Jp, the ITE can be 
simplified to a Poisson equation with Neumann boundary 
condition. 


lb-h ~ V?W(r) + LOWI 


(1) 
Azh+h On 


where fis the focal length, / is the defocus length, Az = oe 


is the conjugate distance of J, r is the coordinate at the pupil 
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Figure 2. Reconstruction of the wavefront was introduced in the Mephisto. First row: intra-focus image, extra-focus image, curvature signal. Second row: curvature 


signal after 5, 10, 20 iterations. 


Table 2 
Examples of Wavefront Reconstruction from Simulated Data 


Zernike Coefficients(unit: A) 


Zernike Term 
Ground Truth1 Estimated 1 Ground Truth2 Estimated 2 Ground Truth3 Estimated 3 Ground Truth4 Estimated 4 

ZA 0.574 0.563 0.204 0.204 0.204 0.196 0.052 0.054 
Z5 —0.601 —0.631 —0.311 —0.293 —0.086 —0.082 —0.327 —0.328 
Z6 —1.009 —1.027 0.222 0.194 0.226 0.216 0.508 0.486 
Z7 —0.680 —0.579 0.830 0.790 0.676 0.651 0.112 0.118 
Z8 —0.172 —0.142 0.204 0.175 —0.649 —0.678 0.512 0.468 
Z9 0.000 —0.013 0.000 0.009 0.000 —0.033 0.000 —0.006 
Z10 0.000 0.014 0.000 —0.017 0.000 0.038 0.000 0.006 
Zll —0.041 —0.039 —0.038 —0.038 —0.038 —0.038 —0.037 —0.038 
Z12 0.001 0.002 —0.001 0.003 —0.002 —0.003 —0.003 0.006 
Z13 0.002 —0.006 0.001 0.001 0.001 —0.002 0.001 0.001 
Z14 0.000 0.001 0.000 —0.009 0.000 0.003 0.000 —0.003 


plane, 7 (r) is the illumination distribution inside the pupil, J, 
and J are the intensity of intra- and extra- focal images, wr) 
is the wavefront, V? is the Laplace operator, 7i is the radial 
direction of the pupil. 

The differentiation of the intensity distribution along the 
optical axis is directly proportional to the wavefront curvature; 
therefore, this technology is commonly called a CWFS 
(Roddier et al. 1988). Various algorithms can be employed to 


estimate the wavefront in curvature wavefront sensing, 
including the Fast Fourier Transform algorithm (FFT) (Roddier 
et al. 1990), the series expansion method, deep learning 
(Patterson & Sutherland 2003), etc. For the popular FFT 
algorithms, the Neumann boundary condition at the pupil edge 
needs to be absorbed into the Laplace operator by setting the 
outside intensity equal to the inside intensity along the radial 


direction n, then ae = 0. Apply FFT transform to both sides 
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of Equation (1), 


1 b-h 27.2 2 = 
FE. = —4r? (u2 + - FLW : 2 


The wavefront can be estimated by, 


1l h-hh 
F{ ze) 


—4r? (u? + v?) 


W(r) =F"! (3) 


As emphasized by Roddier (Roddier & Roddier 1993), the 
Poisson equation serves as a first-order approximation applic- 
able primarily to small Az values corresponding to highly 
defocused images. Furthermore, aberrations must remain 
sufficiently small to confine the distortion of defocused images 
within the linear range (Roddier & Roddier 1993). However, in 
the presence of substantial aberrations, it becomes imperative 
to enhance accuracy through iterative compensation aimed at 
compensating the distortion inherent in defocused images. 
Assuming the system is paraxial, the coordinates (x, y) in the 
defocus plane is a scaled-down of the normalized coordinates 
(xp/R, Y/R) in the pupil plane, that is x= Ix,/Rf and 
y=ly,/Rf, where R is the radius of the pupil. An aberration 
W (xp, y,) in the pupil plane will produce coordinate deviation 
of the rays, 


v =r cy 
x 
(4) 
OW (x, 
y=y+ cW) 
dy 
where C = ees . The intensity of point (x’, y’) and point 
(x, y) has to satisfy flux conservation, 
ow ow 
I(x, y) = (x, y) 1 + C| —— + —— 
(x, y) = 1" nl & a 
2w 82 2w \? 
C2 o W O-W oW l (5) 
Ox? Oy? Oxy 


The widely employed FFT algorithm, coupled with intensity 
compensation, has found extensive application in wavefront 
estimation. In this study, we applied this algorithm to the 
wavefront reconstruction of Mephisto, subjecting it to 
extensive testing using simulated defocused images character- 
ized by known low-order aberrations. Figure 2 depicts intra- 
and extra-focus images simulated through ray tracing of the 
Geometric Image Analysis in ZEMAX. The curvature signal 
exhibited diminishing magnitudes as the defocused image 
intensities were compensated iteratively by the estimated 
aberrations. After 20 iterations, discernible signals became 
negligible. Across all cases outlined in Table 2, wavefront 
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reconstruction consistently yielded errors below 5%. The 
efficacy of the algorithm was further evaluated on the real 
Mephisto telescope, with performance details provided in 
Section 4.2. 


3. Implementation of Multi-field Split Curvature 
Wavefront Sensing 


3.1. The Defocus Distance 


The defocus distance emerges as a pivotal parameter in 
configuring multi-field split Curvature wavefront Sensing. With 
an increase in defocus distance, the doughnut-shaped images 
occupy a larger pixel area, concurrently diminishing the signal- 
to-noise ratio (S/N) and elevating the likelihood of blending in 
sky regions with dense stellar objects. Additionally, the 
intensity patterns of defocused images become less sensitive 
to the wavefront, a consequence of reduced Fisher information 
(Fienup et al. 1998). The advantages conducive to curvature 
wavefront sensing include the weakening of atmospheric 
turbulence effects on the doughnut-shaped images and an 
enhancement in the adequacy of pupil sampling. 

Recognition of the wavefront error through the algorithm is 
contingent upon distorted coordinates of the ray deviation 
exceeding 1 pixel, that is, 


Jara le iar =at e lpixel. 
Ox Oy 


(6) 


Expand the wavefront error in terms of Zernike polynomials, 


Ow OZ; OW OZ; 
= J ZC; d: == > zG ——. 7 
Ox “e Ox a Oy a Oy 


Here zc; is the Zernike coefficient, Z; is theith Zernike 
polynomial. By analyzing the partial derivatives expression of 


ƏZ: _ 
A Ip = r} and 


the Zernike polynomials in Table 3, max { | 


{ | OZ; 
max 


ay 


at different pupil position r can be 


Ip = r} 
calculated (shown in Figure 3). To enable sufficient data for the 
wavefront reconstruction, we take the number of recognizable 
pixels to be greater than half and take the pupil radius of 0.4 as 
a reference. Assuming a wavefront reconstruction accuracy of 
0.01A, then Z4(p = 0.4) ~ 0.48, the upper limit of defocus 
distance should be, 


a(£)y - 0.01 - 0.48 


Iz (8) 


1 pixel 


On the other hand, the defocus images of CWFS must stay 
away from the caustic zone to avoid the intensity patterns 
smeared by the blurred point-spread function. For the 
observatory site to have Kolmogorov spectrum with Fried 


Research in Astronomy and Astrophysics, 24:085006 (15pp), 2024 August 


4 


3.5 


3 


|dZ/dx| 
N 


Ox 


Figure 3. (a) max { 


Table 3 

The Partial Derivatives Expression of the Zernike Polynomials 
Polynomials 
Number Expression Zy Zy 
Z4 (a +y?) —1 4x 4y 
Z5 2y 2y 2x 
Z6 ry 2x —2y 
ZI [3(x? + y?) — 2] 6xy 3x7 + Oy? — 2 
Z8 xB? + y2) — 2] 9x7 + 3y? -2 6xy 
Z9 yx? — y3) 6xy 3(x2 — y3) 
Z10 x(x? — 3y?) 3(x? — y?) —6xy 


parameter ro, the defocused distance / must satisfy (Rod- 
dier 1988): 


= <l (9) 


For example, for the system with parameters of D = 2 m, 
f=4m, A\=0.5 um, lpixel=15 um, ro = 10 cm, the range of 
the defocus distance would be, 


0.2 mm </< 4mm. (10) 


3.2. Signal Pre-processing and Deblending 


In the context of simulated defocused images depicted in 
Figure 4(a), multi-field split curvature wavefront sensors 
capture a bunch of doughnuts simultaneously. To isolate 
optimal doughnuts, the initial step involves binarising the 
intensity pattern of these doughnuts, as represented by the false 
color in Figure 4(a), through the application of a threshold. The 


} at different pupil size; (b)max { 
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4 
3.5 
3 


2.5 


|dZ/dx| 
N 


A 


X 0.40201 
Y 0.482582 


2a } at different pupil size. 


i 
oy 


threshold value is determined using the OTSU algorithm. 
Subsequently, the doughnuts in the binarised image are 
subjected to segmentation and labeling based on connected 
area analysis. Here, the connected area pertains to a region in 
the image characterized by identical values in a 4- or 
8-neighborhood of the pixel. Figure 4(c) presents the 
segmented doughnuts with different colors denoting distinct 
tagged values. By utilizing this segmentation approach, 
doughnuts with a specific S/N can be readily identified and 
selected. 

Due to the expansion of the point source into doughnuts 
consisting of hundreds of pixels, there is an increased 
probability of blending between two adjacent doughnuts in 
densely sky regions. Figure 5(a) illustrates this phenomenon, 
where doughnut images from fainter stars overlay the target 
doughnut, introducing background noise. To mitigate the 
impact of fainter doughnuts on the curvature signal, it is 
imperative to eliminate the intensity contribution of fainter 
doughnuts in the overlapping regions. To achieve this, we 
initially employ the Sobel operator to detect edges in the 
blended image. Subsequently, the blended image is segmen- 
ted into distinct regions based on these edges. By labeling 
these segmented regions, we can accurately identify and 
extract the intensity of the fainter doughnut in non-over- 
lapping areas, fitting the data to polynomials as depicted in 
Figure 5(c). In areas of overlap, the background noise 
originating from the fainter doughnut can be estimated and 
removed through the application of the fitted intensity 
distribution. Figure 5(d) presents the deblended target 
doughnut, with iterative processes contributing to enhanced 
deblending accuracy. 
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3.3. Field-dependent Corrections 


For the split CWFS, the curvature signal is derived from 
intra- and extra-focal doughnuts originating from different 
stars, each possessing non-identical optical properties. 

In these optical properties, the intensity difference can be 
effectively overcome by normalization. Furthermore, the intra- 
and extra-focal doughnuts, sourced from different regions of 
CWES, inherently exhibit non-common path aberrations and 
variations in pupil geometry due to distinct beam footprints 
(Xin et al. 2015). Compensation is applied to mitigate the 
impact of non-common path aberrations by addressing the 
intrinsic design aberrations. Figure 6 illustrates the Zernike 
coefficients representing these intrinsic design aberrations 
across the curvature wavefront sensor field for Mephisto. By 
using the center field of the curvature wavefront sensor as a 
reference, non-common path aberrations can be subtracted 
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Figure 4. (a) Image of the simulated defocus CCD; (b) the binarised image (c)the labeled doughnuts of the image; (d) the picked out doughnut. 


through intensity compensation. This results in the migration of 
doughnuts from different fields to the center field of CWFS. 


3.4. Coordinate Correction for the Intensity 
Compensation 


As described in Sections 2.2 and 3.3, both the iterative FFT 
and the series expansion algorithm for solving the Poisson 
equation with Neumann boundary conditions represent first- 
order solutions. Enhancing estimation accuracy is achievable 
through the integration of the intensity compensation algo- 
rithm. This algorithm assumes a pivotal role not only in 
ameliorating the inherent limitations of the first-order solutions 
but also in compensating for non-identical optical properties 
pertinent to Field-Dependent Corrections. Our previous work 
has expounded upon the correction of nonlinear coordinate 
distortion arising from fast-beam optical properties and off-axis 
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Figure 5. (a) The blended doughnuts of two neighborhood signals; (b) the segmented regions; (c) the polynomial fit of the non-overlapping area; (d) the deblended 


target doughnut. 


effects (Wu et al. 2015). Figure 7 describes the coordinate 
distortion of fast-beam and off-axis mode. 


OW (x,y) 
Ox 

OW (x,y) (11) 
dy 


x = D(x, y) + C 
y = Do(x, y) + C 


D Se eee 
1x, y) VP —@2 +y>R2 x 


D f 
Dit, ) = eae 
E f . : - 
D(x, y) = WP (499K? cos 3— sin pcos 8 
f 1 


z ; R 
f? — (x? + y?)-R2 cosa — sin acos a 


where for the on-axis fast-beam 


system and 


Dy, y) = 


for the off-axis fast-beam system. 


Bo Xin et al. (Xin et al. 2015) introduced a more 
comprehensive approach through numerical mapping, transi- 
tioning from the pupil (x — y plane) to the intra- and extra-focal 
images (x’ — y’ plane). In this method, a set of 10th-order two- 
dimensional polynomials accurately delineates the relationship 
between the two coordinate sets. This generalized numerical 
mapping provides a more precise depiction of the nonlinear and 
non-axisymmetric projection inherent to a specific optical 
system. Using Mephisto as an illustration, we calculate the 
distorted coordinates for the fast-beam mode and off-axis mode 
using both analytical and numerical methods. Figure 8 
illustrates the comparability of distorted coordinates derived 
from these two methodologies, indicating a high degree of 
similarity. 
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Figure 6. Zernike coefficients distribution of the non-common path aberrations. 


P~ 


(a) 


Figure 7. The coordinate distortion of (a) fast-beam mode and (b) off-axis mode. 


4. Multi-field Split CWFS for Mephisto 
4.1. Optical Configuration of Mephisto 


Multi-channel Photometric Survey Telescope (Mephisto) is a 
large field ground-based telescope with a 1.6 m primary mirror 


and 2° field of view (Liu 2019). Figure 9 describes the optical 
layout of Mephisto, three films coated cubic prisms were 
adopted for delivering a 3-channel colored movie of the 
universe. This innovative approach enables the telescope to 
provide real-time, high-precision colors for celestial objects, 
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Figure 8. Nonlinear coordinate mapping of the pupil to defocus plane. (a) Analytic mapping; (b) numerical mapping; (c) the difference between the two methods. 
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Figure 9. Optical layout of Mephisto. 


fostering rapid and robust classification of variables and 
transients. Mephisto’s capabilities extend to delivering panora- 
mic and panchromatic documentation of our evolving universe, 
marking a significant advancement in observational astronomy. 

Driven by considerations encompassing image quality 
requirements, light obscuration, and focal instruments assign- 
ment, the Ritchey-Chretien (RC) optical system, augmented 
with a lens corrector, has been chosen for Mephisto. The 0.7 m 
secondary mirror of Mephisto proves to be notably sensitive to 
alignment adjustments, owing to the fast focal ratio of the 
primary mirror. The field-dependent image quality of mis- 
aligned Mephisto optics manifests as asymmetric, continuously 
varying, and nonlinear with the misalignments. To address this 
challenge, a 5-degree-of-freedom (DOF) mechanism has been 
designed to actively adjust M2 mirror positions based on 


feedback from the multi-field wavefront sensor. This adaptive 
mechanism ensures optimal image quality under operational 
conditions. It is noteworthy that the compact configuration of 
Mephisto allows for a mere 15mm distance between the focal 
plane and the last optical surface. In such a crowded space, the 
CWFES stands out for its ease of installation in comparison to 
the SHWFS, especially in crowded spaces. 


4.2. First Light for CWFS on Mephisto 


Mephisto is being built by NIAOT and installed at Lijiang 
Gaomeigu Observatory (100°2' E, 26° 42’ N, 3200m). On 
2023 June 16, experimental tests of the curvature wavefront 
sensor were conducted. System aberrations were diligently 
minimized, and the camera was positioned at the plane 
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Figure 10. Linearly increase of the introduced coma term in the Mephisto. From left to right: intra-focus image, extra-focus image, curvature signal; From top to 


bottom: decent x = 0 mm, decent x = 0.08 mm, decent x = 0.2 mm. 


corresponding to the best focal plane with minimum stellar 
images. Because the signal area with about 50 x 50 pixels has 
relatively better wavefront sensing results, the intra- and extra- 
defocus images were systematically captured by displacing the 
secondary mirror 0.28 mm along the optical axis, effectively 
simulating a defocused distance of 3.4mm. Employing ray 
tracing of Mephisto, a linear increase in the decentering of the 
secondary mirror introduced an indicated coma value. The 
series of defocus images, illustrated in Figure 10, portray the 
evolving impact of increased decentering of the secondary 
mirror, progressing from 0mm to 0.08 mm and 0.2 mm. The 
first two columns present the intra-focus and extra-focus 
images, and the third column depicts the curvature signal (S). 
Evidently, the intensity accumulates toward the bottom of the 
defocus images with an increase in the decenter x of the 
secondary mirror, concurrently indicating an increase in 
the coma. 


10 


The wavefront is estimated collectively by the inner iteration 
of the Poisson equation solver and the outer iteration of 
intensity compensation. To assess the effectiveness of curva- 
ture wavefront sensing, an examination of the linearity between 
coma coefficients and the decentering of the secondary mirror 
is conducted. Figure 11 illustrates the variation in estimated 
Zernike coefficients corresponding to the increasing decenter- 
ing of the secondary mirror along the x-direction. The 
coefficient of determination serves as a quantitative measure 
of linearity, 

n 42 
R=1— Dini) i, (12) 

ZiO; — Yi) 
Where y; is the Z7 coefficients, J, is the fitted value, y, is the 
mean value of y, n is the sample size. The coefficient of 
determination (R) for Z7 attains a high value of 0.993, 
indicating a strong correlation and affirming the robust linearity 
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0.4 
03 
0.2 
0.1 
x 0 
-0.1 
-0.2 
-0.3 
-0.4 
-0.5 
-0.6 
Z4 Z5 Z6 Z7 Z8 Z9 Z10 ZLE 
—e— dx=0mm 0.2133 -0.0234 0.0914 -0.0584 0.0596 -0.026 -0.0042 0.0459 
—@—dx=0.02mm 0.2336 -0.0678 0.0399 -0.0792 0.0724 -0.0832 -0.0088 0.0402 
—e—dx=0.04mm 0.2193 0.0175 -0.0043 -0.172 0.0633 -0.053 -0.0451 0.0384 
—®—dx=0.08mm 0.2155 -0.0888 -0.0325 -0.2395 0.0549 -0.0299 0.0276 0.038 
—edx=0.12mm 0.2601 -0.0135 0.0035 -0.3583 0.0537 -0.0017 0.0088 0.0342 
—edx=0.l6mm 0.2191 -0.0037 0.0193 -0.4361 0.0681 0.0458 0.0651 0.0404 
—edx=0.20mm 0.2605 -0.1824 -0.0232 -0.5573 0.0604 -0.0446 -0.0538 0.0394 
Figure 11. The estimated Zernike coefficients with different decenter x. 
Table 4 
The Partial Derivatives Expression of the Zernike Polynomials 
Misalignment Parameters 
Iteration M2 M2 M2 M2 M2 CCD CCD CCD 
Piston (mm) Decenter x (mm) Decenter y (mm) Tilt x (°) Tilt y (°) piston (mm) Tilt x (°) Tilt y (°) 
Introduced 
2 2 0 0.02 0.02 0 0.02 0.02 
Estimate 
1 1.27 x 107° 1.74 x 107° —2.47 x 107° -142x104  1.52x 107 = 3.43 x 107? - 1.91 x 1074 3.04 x 1074 
2 2.34 x 1077 6.89 x 1077 7.17 x 1077 —8.40 x 108 83x107% 6.37 x 107’ 2.65 x 1077 —2.50 x 1077 


between the Z7 coefficients and the decentering of the 
secondary mirror along the x-direction. Notably, the Z7 term 
demonstrates linear enlargement, while the Z8 term consis- 
tently remains small, confirming their decoupled nature. 
Consequently, decentering in the x-direction does not influence 
the Z8 term. The Z4 term stabilizes within the range of 
0.21-0.26, indicating a slight offset of the telescope focus from 
the optimal position. Additional confirmation of the linearity 
between Z7 coefficients and decentering in the x-direction is 
presented in Figure 12(a). The difference between the 
simulated and experimental results, as depicted by the error 
bars measuring less than 0.1. The oscillation observed in the 
astigmatism terms Z5 and Z6 can be attributed to the slight tilt 


11 


of the secondary mirror. Figure 12(b) presents an additional set 
of measurements involving a decentered secondary mirror in 
both the x and y directions. Notably, a discernible linearity is 
observed in the Z7 terms concerning decentering in the x- 
direction and the Z8 term in relation to decentering in the y- 
direction. Figure 13 is a comparison between the theoretical 
adjustment amount of the secondary mirror and the actual 
measurement value. 

We introduced an unknown decentering to the secondary 
mirror, subsequently measuring the coma aberration through 
the curvature wavefront sensor. With the linear relationship 
depicted in Figure 12, the decentering value of the secondary 
mirror is determined and corrected. As shown in Figure 14, 
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Figure 12. (a) The linearity of Z7 as the decenter x increases; (b) the linearity of Z7 and Z8 as decenter x and decenter y increases. 
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Figure 13. Comparison of theoretical adjustment and the actual measurement. 


even with large fluctuations of the seeing and incomplete 
calibration of other misalignments (FWHM about 2”), the 
images captured after introducing the decentering of the 
secondary mirror exhibit noticeable coma aberration. After 
correction, the coma aberration has been effectively mitigated. 


4.3. Performance of Active Alignment for Mephisto 


An application has been specifically developed for Mephisto, 
with the user interface (UI) showcased in Figure 15. Within this 
software, four fundamental modules are defined: 


1. CWFS Control: This module facilitates the configuration 
of parameters related to the CWFS, encompassing optical 
characteristics of Mephisto (diameter, focal length, 
obscuration), CCD parameters (pixel size, pixel number), 
and defocus distance. 


2. Source Selector: This module enables the selection of 


source signals for the wavefront sensors and calculates 
the field associated with each signal. 


. Wavefront Estimator: Source signals are transposed from 


their original field positions to the central field position of 
the CWFS. Subsequently, the wavefront is estimated 
through the inner iteration of the Poisson equation solver 
and the outer iteration of compensating for estimated 
aberrations. 


4. Misalignment Calculator: This module determines the 


rigid body motions of the second mirror and the focal 
plane through a sensitivity matrix. 


The fundamental functions of the first three modules have 
been illustrated in Section 3. Here is a comprehensive 


description of Module 4. The relationship between misalign- 


ments and wavefront errors of each field point can be 
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(b) 


Figure 14. (a) The image captured after introducing the decenter of the secondary mirror; (b) the image captured after correction of the decenter of the secondary 


mirror. 


represented by, 


Fi A, X2 Xn) Foi 
2 = h X2, Xn) ca ie (13) 
Fin Sin xı, X2,°° Xn) Fom 


Here F; is the Zernike coefficients of wavefront errors at each 
field point, Faj is the baseline design aberrations, function fj 
characterizes the asymmetric, continuous varying and nonlinear 
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relationship, x; is the misalignment parameters. Finding an 


exact expression for f; proves to be challenging; nevertheless, f; 
can be effectively approximated through linear equations, 
particularly when subjected to minor misalignments. 


Of; j 
F; = Fo + — a — xo) + —@2 — x02) 
Ox, 0. 


X2 
; (14) 

Of; 

se des +5, C” — Xon) 


Xn 
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Figure 15. The UI of the active optics for Mephisto. 


where 


a is the partial derivative of the Zernike aberration f; to 


the misalignment parameter x; which can be approximated as 
difference equations, 


7 of 
F; = Foj + Loa — X01) + — (x — X02) 
bx, 6x2 


of. 
Spe die = Xon). (15) 


OXn 


Then the approximate linear equations for aberration and 
misalignment are, 


KH Hi of 
F F Ox, 6x2 OXn 
nales. al poa 
: 6x, 6x2 OXn : . (16) 
Fn = Fom Xn — Xon 
fn Sfm fn 
6x, 6x2 OXn 
Represent the equations in matrix, 
AF =S. Ax (17) 


where AF is the deviation of the measured aberration from the 
designed aberration, Ax is the value of each misalignment 
parameter, S is the sensitive matrix. 

S is typically computed by numerically introducing pertur- 
bations with specific step lengths (Ax, Ax,---,Ax,) into the 
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optical system. The choice of step length for each perturbation 
is contingent upon its sensitivity, a metric assessable through 
the eigenvalues of S. Perturbations characterized by larger 
eigenvalues, such as tilt-x and tilt-y of the second mirror, 
indicate a heightened sensitivity, necessitating smaller step 
lengths to confine aberrations within the linear range. 
Conversely, less sensitive perturbations, like the tilt-x and 
tilt-y of the camera, require larger step lengths to ensure the 
magnitude of aberration changes surpasses wavefront sensing 
errors. 

The misalignment parameters Ax can be solved from the 
measure wavefront errors, 

Ax = (STS) ISTAF. (18) 

In this study, the strategy employed for simulating active 
alignment involves utilizing Dynamic Data Exchange (DDE) 
between ZEMAX and MATLAB. Perturbation values are 
generated through MATLAB and subsequently transferred to 
ZEMAX. The wavefront errors for each field point are then 
determined through ray tracing and conveyed to MATLAB for 
the calculation of the sensitive matrix and compensator signals. 
The retrieved misalignment parameters from simulated data are 
presented in Table 4. Perturbations corresponding to the rigid 
body motion of the second mirror (M2) and CCD are 
introduced with their respective values outlined in Table 4. 
Importantly, the remaining perturbations after a single iteration 
fall within the tolerances defined for the telescope. 
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5. Conclusions 


This work presents a comprehensive analysis of the multi- 
field split curvature wavefront sensing system. The upper limit 
of the defocus distance is determined based on the distorted 
coordinates of ray deviation. The labeling of doughnuts from 
various sources involves an initial categorization by connected 
areas, facilitating the subsequent selection of optimal signals. 
For deblending overlapped doughnuts from neighboring 
sources within the blended image, different segment regions 
are defined. The non-overlapping areas of fainter doughnuts are 
subjected to polynomial fitting, and overlapping areas are 
estimated and removed through these fitted polynomials. 
Compensation for non-identical optical properties between 
intra- and extra-focal images is implemented. Both analytical 
and numerical methods are employed to calculate distorted 
coordinates for fast-beam and off-axis modes, with a 
subsequent comparative analysis of the results. 

The curvature wavefront sensor was effectively tested on 
Mephisto, demonstrating its capability to diagnose telescope 
aberrations accurately. Specifically, we conducted tests focus- 
ing on the prevalent coma aberrations arising from the 
decentering of the second mirror. The observed coma exhibited 
consistent linearity with the incremental decentering of the 
second mirror. Furthermore, simulations were performed to 
assess the active alignment performance of Mephisto, yielding 
retrieved misalignment parameters from simulated data that fall 
within the defined tolerances for the telescope. 
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